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We derive and systematically analyze scalar glueball correlation functions in both the hard-wall 
and dilaton soft-wall approximations to holographic QCD. The dynamical content of the holographic 
correlators is uncovered by examining their spectral density and by relating them to the operator 
product expansion, a dilatational low-energy theorem and a recently suggested two-dimensional 
power correction associated with the short-distance behavior of the heavy-quark potential. This 
approach provides holographic estimates for the three lowest-dimensional gluon condensates or al- 
ternatively their Wilson coefficients, the two leading moments of the instanton size distribution in 
the QCD vacuum and an effective UV gluon mass. A remarkable complementarity between the 
nonperturbative physics of the hard- and soft-wall correlators emerges, and their ability to describe 
detailed QCD results can be assessed quantitatively. We further provide the first holographic esti- 
mates for the decay constants of the 0"*"'" glueball and its excitations. The hard-wall background 
turns out to encode more of the relevant QCD physics, and its prediction fs — 0.8 — 0.9 GeV for the 
phenomenologically important ground state decay constant agrees inside errors with recent QCD 
sum rule and lattice results. 



I. INTRODUCTION 



Despite more than three decades of intense experimental and theoretical scrutiny the long predicted glueball states 
[l| of Quantum Chromodynamics (QCD) remain stubbornly elusive 0, The slow pace of theoretical progress 
reflects the extraordinary complexity of the infrared Yang-Mills dynamics which generates both the gluonic bound 
states and their mixing with quarkonia. New analytical approaches for dealing with strongly coupled gauge theories, 
as they have recently emerged from gauge/gravity generalizations of the AdS/CFT correspondence 0j should 
therefore find rewarding and much needed applications in the glueball sector. 

Until now such applications have focused on the glueball mass spectra, which were among the first holographically 
calculated observables in a variety of more or less QCD-like gauge theories @ (for a review and current developments 
see e.g. Refs. [7]). More recently, glueball spectra were also obtained in the first bottom- up [tI] proposals for the 
holographic QCD dual [12, [H, Hj) HBl as well as in back-reacted models [26l. [27|. 

In the present paper we are going to extend the holographic analysis of glueball properties beyond the spectrum, 
by focusing on the gauge physics content of the glueball correlation function and its spectral density. We will relate 
the holographic predictions to QCD information from the operator product expansion (OPE), a low-energy theorem 
based on the anomalous Ward Identity for the dilatation current, and a recently advocated, effective UV gluon mass. 
Our calculations will be based on two alternative AdS/QCD backgrounds, namely the AdSs geometry with a "hard 
wall" IR brane cutoff (of Randall-Sundrum type f2^) in the fifth dimension [8] and the dilaton-induced soft wall p^ . 
which both proved phenomenologically successful in the meson sector [9, V2, 23] . 

A second major objective will be to provide the first holographic estimates for the decay constants of the scalar 
glueball and its excitations, i.e. for the glueball-to-vacuum matrix elements of the lowest-dimensional gluonic QCD 
interpolator. These on-shell observables are of particular interest because they contain fundamental information on 
glueball structure and govern the spacial extent of the glueball (Bcthc-Salpeter) wave functions. Lattice indications 
for an exceptionally small size of the lowest-lying scalar glueball ^22], for example, should translate into an unusually 
larg e value of its decay constant. Evidence for such an enhancement was indeed found in instanton vacuum models 
[30j as well as in those QCD sum rule analyses which include instanton contributions to the OPE coefficients [3ll[33|. 

The decay constants, which are the first glueball observables besides the low-lying spectra for which direct 
(quenched) lattice results are now available [33| , also play a crucial role in the theoretical analysis of glueball pro- 
duction and decay rates. For this reason, their accurate prediction will be instrumental in eventually meeting the 
two longstanding challenges of glueball physics, i.e. the establishment of unambiguous glueball signatures and their 
experimental identification. As a case in point, the decay constants provide critical nonperturbative input for the cal- 
culation of glueball production amplitudes in the "gluon-rich" radiative heavy-quarkonium decays which are currently 
measured at EES 34]. 

The paper is structured as follows: in section |IT] we define the dual bulk dynamics on which our study will be based, 
and we derive general expressions for the scalar glueball correlator and the decay constants in IR-deformed AdSs 
duals with a nontrivial dilaton background. In section Hill we focus on the two AdS/QCD backgrounds mentioned 
above (i.e. hard and soft wall) and derive exact analytical expressions for the corresponding correlators and their 
spectral functions. Wc then analyze the results by confronting them with the OPE of the QCD correlator (including 
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nonperturbative contributions to the Wilson coefficients), the dilatational low-energy theorem which governs its low- 
momentum behavior, and the contributions of an effective UV gluon mass. This strategy will allow us to obtain 
holographic estimates for various QCD vacuum scales, i.e. three gluon condensates (or alternatively their Wilson 
coefficients) and the two leading moments of the instanton size distribution, as well as for an effective UV gluon 
mass. In section IIVI we obtain holographic predictions for the (ground and excited state) glueball decay constants 
and compare them to other available theoretical results. Section |Vl finally, contains a summary of our findings and 
presents our conclusions. 

II. DUAL DYNAMICS OF THE SCALAR GLUEBALL 

The gauge/string correpondence [1, Q maps string theories in curved, ten-dimensional spacetimes into gauges 
theories which live on the d dimensional boundaries. For UV-conformal gauge theories like QCD with d — 4, the dual 
spacetime metric factorizes into a five-dimensional non-compact manifold which close to its boundary approaches the 
anti-de Sitter space AdSs (R) of curvature radius R, and a five-dimensional compact Einstein space (where e.g. 

= (R) for the maximally supersymmetric gauge theory) with the same intrinsic size scale. The corresponding 
line element is Q 

ds^ - gMN (x) dx^'dx^ = e'^(") ^ {ri^Jx^'dx" - dz'') + R^ds\^ (1) 

(in conformal Poincare coordinates) where 77^^^ is the four-dimensional Minkowski metric. Conformal invariance of 
the dual gauge theory in the UV requires the absence of AdS deformations (i.e. A{z) ^ 0) as z ^ 0. Supergravity 
solutions suggest the additional presence of a nontrivial dilaton background <I> (a;) , and potentially of other background 
fields (including Ramond-Ramond axions, tachyons etc., see e.g. Ref. [13, HBl) which do, however, not play an explicit 
role in the AdS/QCD duals considered below. 

A. Bulk action and holographic glueball correlator 

The scalar QCD glueballs are interpolated by the lowest-dimensional gluonic operator carrying vacuum quantum 
numbers, 

Os {x) = [x) G^^^^ (x) , (2) 

(where C^^ is the gluon field strength) which also figures prominently in the anomalous dilatational Ward Identity 
and in the corresponding low-energy theorems (cf. Appendix |A| . Since the conformal dimension of Os is A = 4 
(at the classical level), the AdS/CFT dictionary prescribes its dual string modes tp{x,z) to be the normalizable 
solutions of the scalar wave equation in the bulk geometry ([Ij (and potentially other background fields) with the UV 

behavior tp {x, z) z^cf) (x). The latter implies that the square mass m\R^ = A (A — rf) = of the bulk field 
Lp vanishes, and that its minimal action has the form 

S [p- g^^^±^j d^+'x^\e~''g^'''dMVdNV (3) 

(where can be related to the five-dimensional Newton constant [s^) which we decompose for later use as 5 = 
Sm + Sqm into bulk and boundary contributions with 

SMW;gM = -^ I rf'^+^xv^e-*^ [e*VMe- V'^a^] (4) 

(Vm = Om + l^r^^^ Qm IsT^^) and 

SaM 9,^]^7^ f [a^ (z) e-'^pd,^] (5) 

JdM 

where (z) = exp2A (z) is the warp factor. The boundary dM consists of the UV brane z = e ^ and of 

an additional IR brane at z — z„i in the hard-wall geometry. 

Variation of the bulk action ^ with respect to p yields the field equation 

e*VMe-*ff*^^5w^ {x, z) = [A ~ {dM<^) g^'"" On] P [x, z) = (6) 
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where A = Va/V is the Laplace-Beltrami operator of the metric ([T]). The action density of the solutions is finite 
only on the boundary dM while S'|°"~'''^°'''' = 0. We now specialize to dilaton fields <I> which depend exclusively on 
the fifth dimension, i.e. = (z). The d-dimensional Fourier transform (p (g, z) of the normalizable dual modes then 
solves the reduced field equation 

[dl + (d - 1) {a-^d,a) ~ + q^] (p {q, z) = (7) 

with discrete on-shell momenta [75] — rn^ in both hard- and soft-wall backgrounds. The eigenvalues m'^ determine 
the glueball mass spectrum of the boundary gauge theory, and the corresponding orthonormalized solutions will be 
denoted ip„ (z) = N^^p {rrin, z). 

Holographic glueball correlation functions are obtained by differentiating the bulk action of the solutions with respect 
to the boundary source . The on-shell action can be constructed with the help of the bulk-to-boundary propagator 
K {q, z) 0, which is the solution of the field equation ([7]) subject to the UV boundary condition K [q^e ^ Q) = 1. Its 
spectral representation is therefore 

n 

(where the limit e — > at the end of the calculation is implied) and automatically satisfies the IR boundary condition 
imposed on the tp^ (z). Hence one can write the solution of Eq. ([6]) corresponding to a given boundary source tp^'^'^ [x') 
as 

^ (x, z) = / T^e-'^-X {q, z) f d^x'e^^-'^^^^ {x') (9) 
J (27r) J 

and obtain the associated on-shell action (which plays the role of a generating functional) by inserting Eq. ([9]) into 
the surface action Taking two functional derivatives with respect to p'-"^ then yields the two-point correlation 
function 



{TOs (x) Os (y)) = I I -^e-'^(--^)n (10) 
(27r) 



of the scalar glueball where 
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-Kiq,z)d,Kiq,z) 



(11) 



For $ fields which vanish at the UV boundary (as does the soft wall dilaton considered below), furthermore, the 
bulk-to-boundary propagator in the form ([5]) generates the spectral representation 

f^r „2N y y (e) V'n (e) _ y fnK .^2) 

^ ' \ KE^ J ^ q^ — ml -I- ie ^ g2 — -H is 

of the correlator (where a prime denotes differentiation with respect to z and (divergent) contact terms are not written 
explicitly). The pole residues of Eq. at = ml define the decay constants of the n-th 0++ glueball excitation 
as 

/,.:=J,(0|O.(0)|0,r> = ^*ya. (13) 

The physical role of the decay constants as the glueball "wave functions at the origin" becomes more transparent 
when viewing them as the coincidence limit of the Bethe-Salpeter amplitudes 

X„(.)^(0|2tr{G,.(-|)^(-|,|)G-(|)}|0-) (14) 

(where the adjoint color parallel transporter U {x,y) ensures gauge invariance and proper renormalization of the 
operators is understood). A smaller glueball size implies a higher concentration of the wave function at the origin 
and consequently a larger value of /„. Since the decay constants are on-shell observables related to the bilinear part 
of the bulk action, we expect them to be reasonably well predicted by Eq. even though the dual dynamics ^ 
contains operators of minimal dimension only. 
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B. Comments on the scalar dual dynamics 

We have restricted the action ([3]), i.e. the dynamics of fluctuations dual to scalar glueballs in a metric plus dilaton 
background, to contain only operators with the minimal number of fields and derivatives. This is appropriate for the 
AdS/QCD candidates under consideration and will permit us to derive analytical expressions for the exact correlators. 
These restrictions also entail several typical limitations of contemporary bottom-up models, however, which we now 
discuss in view of their potential impact on the glueball sector. 

A first obvious limitation is the treatment of the metric as a non-dynamical background field, reflected in the fact 
that the action ^ contains neither the Einstein-Hilbert term nor higher-derivative corrections to it [t^. Graviton 
fluctuations around the bulk metric are generally neglected as well in bottom-up duals. If included, they could have 
a direct bearing on the scalar glueball dynamics since fluctuations around a non-conformal metric would generate a 
scalar "radion" mode (related to the fifth or radial dimension component of the graviton) as it appears e.g. in distance 
fluctuations between the branes of Randall-Sundrum models ^28l, Isoj and in dynamical dilaton-gravity models (see 
e.g. Ref. ^26J). By mixing with the scalar held in the bulk action ([3]) the radion would then modify the glueball mass 
spectrum and the diagonalized correlator. Since the mixing strength decreases with increasing mass of the excitation 
dual to the radion (or mixed radion-scalar), the standard neglect of the graviton dynamics e.g. in hard- wall models 
is equivalent to the tacit assumption that an a priori unspecifled stabilization mechanism pushes this mass up far 
enough for radion admixtures to become negligible (ideally by breaking conformal symmetry according to the QCD 
trace anomaly). 

A more realistic holographic description of the glueball sector would probably require to include further operators, 
potentially of higher dimension, into the action ([3]). Those may contain additional bulk flelds, with promising candi- 
dates including spin-zero background fields encoding condensates of relevant QCD operators and flavor-carrying gauge 
flelds for the description of quarkonium-gluonium mixing effects and for speciflc glueball decay channels (see e.g. 
Ref. [37|). Operators containing a higher number of derivatives are potentially important as well. They typically arise 
from stringy a' corrections in bulk regions where the curvature radius R of the geometry becomes comparable to the 
string length Ig . In holographic duals of large- iV^ gauge theories such regions are expected to describe the UV regime 
where the 't Hooft coupling A — gyj^Nc becomes small, i.e. where A ~ [R'^/a') < 1 [t^]. The lack of asymptotic 
freedom in current bottom-up duals (as well as in supergravity approximations), i.e. the fact that AdS/QCD models 
remain strongly coupled in the UV (although they approach a conformal flxed point) 78], is therefore closely related 
to the absence of higher-dimensional operators. 

This discussion indicates that improvements of the AdS/QCD approach will depend in no small measure on whether 
a quantitative understanding for the impact of the strongly coupled UV regime on holographic predictions can be 
developed. In the present paper we propose a strategy towards clarifying this issue which is based on the comparison of 
holographic model predictions for hadronic correlation functions with the QCD operator product expansion (OPE). 
The OPE lends itself particularly well to a systematic diagnosis of the UV sector since it factorizes gauge theory 
amplitudes into short-distance mode contributions to the Wilson coefficients and long-distance physics contributions 
to local operators. We are going to exploit this factorization property below when searching for speciflc traces of the 
strongly coupled UV regime in the two-point function of the scalar glueball channel and will indeed find evidence 
for deficiencies in the AdS/QCD description of the perturbative Wilson coefficients (beyond the leading conformal 
logarithm). Moreover, the results will suggest systematic improvement strategies for bottom-up duals. 

In view of the issues raised above one may wonder whether a decent holographic description of asymptotically free 
Yang-Mills theories at large Nc could at all be achieved on the basis of a local flve-dimcnsional action (which may 
include a few higher-dimensional operators). Fortunately, there are several indications for an afflrmative answer. A 
general argument due to Witten implies that the locality of the five-dimensional bulk dynamics is ensured by the large- 
Nc limit [4^ . The bulk action may then be viewed as an effective string field theory which contains an elementary field 
for each string excitation (including those of arbitrarily high spin) while higher-dimensional operators are suppressed 
by powers of 1/Nc [l^]- Moreover, a' corrections may be partially resummed e.g. into the dilaton potential [27j . 
and extensive QCD sum-rule [1^ analyses have shown that already a few leading OPE power corrections, and hence 
hopefully the few corresponding light bulk fields with controUably small a' corrections, can capture at least the 
essential properties of most hadronic ground states. 

III. HOLOGRAPHIC GLUEBALL CORRELATORS 

The expressions derived above hold for all geometries of the form ^ and for general dilaton backgrounds $ {z) 
with $ (0) — 0. In order to gain dynamical insight into the holographic glueball correlator and to obtain quantitative 
estimates for the decay constants, we will now consider two s pec iflc AdS/QCD backgrounds, i.e. the AdS^ slice of the 
hard IR wall geometry and the dilaton soft wall of Ref. jl2l|. In particular, we will derive analytical expressions 
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for the glueball correlator ((TT|) and its spectral density in the hard- and soft-wall backgrounds. Those will then be 
analyzed by comparison with the QCD operator product expansion, a dilatational low-energy theorem which governs 
the correlator at zero momentum, and an effective UV gluon mass contribution of the type suggested in Ref. [31 ■ 
The pertinent QCD information is summarized in Appendix [Xl 



A. Conformal symmetry breaking by an IR brane 



A substantial part of the successful recent AdS/QCD phenomenology (see e.g. [8|, |9|, [l^, |23| ) was obtained on the 
basis of the so-called "hard wall" geometry \§\ . This rather minimal deformation of the AdS5 metric approximately 
describes IR effects including confinement by a sudden onset of conformal symmetry breaking in the form of an IR 
brane at z = Zm, i.e. 

e'^''"'(^) =0(z„-z), zm^^^hu^ ^C^-^^O, (15) 

which reduces the five-dimensional bulk spacetime to an AdSs slice. 

We start our analysis of the holographic glueball correlator in this highly symmetric background where analytical 
expressions are straightforward to obtain. The bulk-to-boundary propagator K [q, z), in particular, can be found by 
solving the field equation ([7]) in the geometry (fT5l) . subject to the UV boundary condition K {q;e) — 1 (with e — 0) 
and the Neumann IR boundary condition d^K {q; z„j) ~ 0. The result is [7§| 



K{q,z) = ^ {qz f 



— -J2 {qz) - Y2 (qz) 

•h yqzm) 



(16) 



where Ju,Yi, are Bessel functions and K{Q,z) ~ 1. After plugging Eq. (fT6|) into the general expression (fTTjl and 
analytically continuing to spacelike momenta ~ ^9^, one ends up with the hard- wall glueball correlator 



KliQZra) , (Q 



2 :7" '7 - In, , 
h[QZm) \M 



t2 



(17) 



{Ky,Ii, are McDonald functions [45[) where two contact terms associated with UV divergent subtraction constants 
were discarded. 

It is instructive to find the spectral density p (s) of the correlator (fTT]) which is defined by means of the dispersion 
relation 

ii{Q^)^rds^ (18) 

(where the necessary subtraction terms are again implied but not written explicitly) and can be derived, e.g., from 
the well-known analyticity properties of the McDonald functions [ist (and the causal pole definition) as the imaginary 
part of II/tt at timelike momenta. The result is 

p(^)™^^f:^fe^ (19) 

where we encountered the hard-wall mass spectrum to„ = ji^n/zm (cf. Eq. P^ ). The spectral weight (jTH]) is non- 
negative, in agreement with general principles, and consists of a sum of zero-width poles, as expected at large Nc 
where the (infinitely many) glueballs become stable against strong decay. The leading large-s behavior of the density 
(fT9|) necessitates subtractions in Eq. ([T8|) and ensures the leading logarithmic dependence of the correlator (fT7|l . 

The holographic result (|17p can be compared to the QCD short-distance expansion (|Aip for Q ^ /i^ > z"^. A 
standard procedure @ for fixing the overall normalization B? /k^, in particular, is to match the coefficients of the 
leading conformal logarithm in Eqs. (|A2p and (fT7|) . which yields 



i?3 _ 2{N'l 

*r2 /TT- 



(20) 



(For a discussion of the accuracy of such estimates see Ref. [1^.) Below we will specialize Eq. (|20)) to the phenomeno- 
logically relevant Nc = 2> which seems - at least as far as glueball properties are concerned - to be a surprisingly 
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good approximation to large Nc [13]. We now turn to the remaining term in the holographic correlator (fT7|) . which 
describes nonperturbative contributions of the boundary gauge theory and becomes 

^ ' 4k2 /i (Qz„) 



TT 



l4g-2Qz„ (21) 



in the OPE limit ^ ^qcd ^m^- ^1^- IpTj) reveals that the hard- wall glueball correlator contains no power 
corrections and that all of its nonperturbative content has an exponential dependence (times powers of Q^). In 
the OPE (|A1|) this exponential behavior originates from small-size instanton contributions to the Wilson coefficients. 
Indeed, for Q ^ the direct instanton contribution (IA6P becomes 



which has exactly the momentum dependence of the first subleading term in the non-perturbative hard-wall correlator 
(I^Tjl . As shown in instanton vacuum models [s^l and directly from the lOPE in QCD sum rules [s!], [s^l, these 
instanton-induced correlations are attractive and of relatively short range ~ p. Hence they reduce the mass and size 
of the scalar glueball while increasing its decay constant. 

For a quantitative comparison of holographic and instanton-induced contributions one may approximately equate 
Eq. ([^ with the second term in Eq. (PTjl . This yields the expressions 

3 1 

(23) 



24527r2C z; 



for the average instanton size p and the overall instanton density n in terms of the hard- wall IR scale Zm- The 
relation p ~ z^, is consistent with the duality between gauge-theory instantons of size p and pointlike bulk objects {D 
instantons or D (—1) branes in the supersymmetric case 48]) localized at a distance z = p from the UV boundary. 
However, it also identifies the instanton's average size p with the maximal size z^ in the AdSs slice, which is likely to 
result in an overestimate. Indeed, the standard identification z"^ ~ Aqcd — 0.33 GeV would imply p ~ 0.6 fm, i.e. 
almost twice the ILM value Pilm ~ 0.33 fm [4§|. As a consequence, niLM — 0.5 fm"'* [4§| would be underestimated 
by the second relation in Eq. 

Besides other likely limitations of the hard- wall background including the strongly coupled UV dynamics, the large 
estimate for p may also reflect the absence of fundamental quark flavors in the simple dual dynamics ([3]). The relations 
(I23p (as well as other results below) may therefore apply more accurately to pure Yang-Mills theory for which several 
lattice studies indeed find larger average instanton sizes p ~ 0.4 — 0.5 fm [50|. In any case, one would not expect 
the instanton scales of the QCD vacuum to be precisely encoded in the hard-wall approximation. In fact, it seems 
remarkable that this minimal background can even semi-quantitatively reproduce the key instanton contribution to 
the short-distance expansion. For a fully quantitative study of such corrections one should resort to top-down gravity 
duals in which the relation between bulk and boundary instantons can be traced exactly [4^. Such investigations 
may also shed light on the interpretation of the leading exponential contribution to Eq. (|2ip in terms of gauge theory 
physics. 

It is interesting to confront the hard- wall correlator with the QCD low-energy theorem (|A7[) . The correlator (fTT)) 
vanishes at = since the removal of the contact terms amounts to subtractions at = 0. Even the contact 
terms do not contain a finite (or infinite) contribution to H(0), however, and neither does the nonperturbative 
part (|2ip alone which would remain after subtracting the perturbative contributions from the spectral density, as 
suggested in the original definition (Slj . (This is in contrast to the one- instanton contribution (jA6[) which contains a 

subtraction term n(^^^) (0) = 2'^5'^(^fi. The one-instanton approximation is not reliable at small Q^, however, where 
multi- instanton and other long- wavelength vacuum field contributions are likely to dominate the correlator.) From 
the LET perspective this is consistent with the absence of power corrections and gluon condensates in the hard-wall 
background. As a consequence, both sides of Eq. (|A7[) vanish identically and the LET is trivially satisfied. A more 
complex situation will be encountered in the soft-wall background below. 

The absence of condensate effects in the hard-wall approximation is not surprising because their purely geometrical 
encoding is known to require power-law deformations [13] of the warp factor A{z) in the infrared [soj . Since large 
instantons generate finite gluon condensates, this furthermore indicates that the correlator ([TT]) receives small-size 
instanton contributions only, in perfect agreement with our above discussion which indeed implies 

P^'^"^ < ^™ ~ p-^ (24) 
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because the instanton size cannot exceed the extension of the AdSs shce in the fifth dimension. Hence the simple 
hard-wall approximation seems to capture the fact that an essential part of the nonperturbative contributions to 
the 0++ glueball correlator is hard compared to the OPE scale fj, and therefore resides in the Wilson coefficients 
[31I [33 , UM- Since the power corrections of the OPE (|A2[) are suppressed by unusually small Wilson coefhcients, 
furthermore, the hard-wall background may indeed provide a reasonable first approximation to the scalar glueball 
correlator. 



B. Dilaton-induced conformal symmetry breaking 

A well-known shortcoming of the hard-wall background (|15p is that it predicts sq uared hadron masses to grow 
quadratically with high radial, spin and orbital excitation quantum numbers fldi llll. [2^. in contrast to the linear 
trajectories expected from semiclassical flux tube models [52]. This problem manifests itself also in the hard- wall 
glueball spectra (cf. Eqs. (|4T|) . (|42)) ) which do not reproduce the expected linear Pomeron trajectory [ssl. [s^]. 

The presence of a nontrivial dilaton background field $ (z) cx was recently proposed as an economical remedy 
for this problem in the meson [12] and glueball |25| sectors. In the simplest version of the resulting gravity dual, 
conformal symmetry breaking in the IR is the exclusive task of the dilaton while the geometry ([T]) remains undeformed 
AdSs, i.e. 

A^'""^ (z) = 0, ^^''^^ (z) = A^z^. (25) 

In the present section we derive and analyze the scalar glueball correlator in this "dilaton soft-wall" back grou nd. 
(Alternative holographic realizations of linear trajectories have been obtained for mesons or glueballs in Refs. [SSIIsgI. 
[53, [H] and for both mesons and baryons in Ref. [Tsl].) 

Although the background (|25p is somewhat more complex than the minimal hard- wall geometry (jlSp . one can still 
find a closed integral representation for the corresponding scalar bulk-to-boundary propagator ([5]) , 



2/2 



l)^'dx(l-a;)x-[«'/(4^')+i]e-T^^ 



Kiq;z)^^l^^-1] I dxil-x) x'i'^ J+^g-T^^ ^ , (26) 

which may be rewritten in terms of confluent hypergeometric functions. Eq. ([26]) is easily shown to be the solution of 
the field equation ([7]) in the background P5|) which satisfies the UV boundary condition K {q; 0) = 1 and additionally 
K (0; z) = 1. Inserting the expression (pS)) into Eq. pT|) leads to 



(27) 



which is the exact soft-wall correlator at spacelike momenta q^ — — Q^. The remaining integral can be performed 
analytically. This is conveniently done by absorbing the small-e singularity into the integrand such that the branch 
cut structure becomes manifest. One then obtains 

ft (0") (S + + i=s + ^- ''^''^ 



where U{a,b,z) is the (multi- valued) confluent hypergeometric function ^45]. After taking the e ^ limit and 
discarding two divergent contact terms, one finally ends up with 



(29) 



in terms of the digamma function (z) ~ T' (z) /P (z) [isj . 

As in the hard-wall case, we begin our analysis of the correlator (|29p by deriving its spectral density from the 
dispersion relation (|18p (where the lower boundary of the integration region is now Smin = ttt-o: see below) as the 
imaginary part of H/tt at timelike momenta. The analyticity structure of the digamma function f45| and the causal 
pole definition then imply 
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The spectral density ((30| is non-negative for s > TOq/^ ^'^d consists, as it's hard-waU counterpart ([T9|) and as expected 
at large N^, of a sum of zero-width poles at the soft-wall masses = 4 (n -I- 2) (cf. Eq. (US])). The leading large-s 
behavior again encodes the conformal large-Q^ behavior of the correlator. 

In order to compare the holographic soft- wall correlator to the OPE (jAip at S> ^qcd^ '^^ rewrite Eq. ([^5]) 
for g2 > 4A^ by means of the asymptotic expansion for the digamma function f45l| and the Bernoulli numbers 
B2n = (—1)"^^ 2 {2n)l( (2n) / (27r)^" {( (z) is Riemann's zeta function) as 



n(Q^) - -^A 



01 

AX- 
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4A2 



2nN 



(31) 



(In the last line we have adapted the correlator to the OPE scale fi by absorbing additional, finite pieces into the 
contact terms). Note that the expansion coefficients grow factorially with the power of \^/Q^, as expected from 
QCD. The coefficients of the conformal logarithm in Eq. (pij) and in the hard-wall correlator ([T7]) are identical. 
This is because large momenta Q probe the 2: — > region where neither the dilaton nor the IR brane affect the 
correlator, so that the same AdSs-induced logarithm governs its behavior in both hard- and soft-wall backgrounds. 
Hence comparison with the perturbative gluon loop of the OPE (jA2ll fixes the normalization / as in Eq. ((20|) 
and as anticipated in the second line of Eq. pTjl . 

In addition to the leading conformal logarithm, the expansion (j31|) contains an infinite tower of power corrections. 
Comparison with the OPE (|A2[) suggests them to be related to the gauge-theory condensates 



{Od) 



(32) 



of = 4,6,8, ... dimensional (local, gauge-invariant) composite operators. The appearance of the scale factor A 
shows that the soft-wall power corrections are entirely dilaton-induced, in contrast to those arising from (hadron 
channel dependent) deformations of the metric in the geometric approach or from additional scalar background 
fields. Tentatively equating the coefficients of the i? = 4, 6 and 8 terms (without O (ag) corrections) to their OPE 
counterparts in Eq. (|A2|) allows for a more quantitative check of the holographic expansion (131]) . The resulting 
relations for the three lowest-dimensional gluon condensates (defined at the OPE scale ~ 1 GcV) are 
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Btt' 



-A« 



(33) 
(34) 
(35) 



These holographic estimates indeed reproduce the order of magnitude expected from QCD. This is mostly because 
their scale is set by the dilaton IR parameter A ^ V2Aqcd uM which generates the mass gap and because the 
coefficients in Eqs. (|33l) and (p4|) are more or less of order unity. The sign of the most reliably determined four- 
dimensional QCD gluon condensate (G^) ~ 0.4 — 1.2 GeV* is positive, however, in contrast to Eq. ([55)1 . QCD 
estimates of both signs exist for the three-gluon condensate, namely the lattice prediction (gG^) ~ —1.5 (agG^)^^^ 
1591 1 and the sing le-instanton value {gG^) ~ 0.27 GeV'^{asG^). The signs of Eqs. ((33| and p5|) . furthermore, are at 
odds with the factorization approximation (60j 



(36) 



for the four-gluon condensate combination (IA4[) . These shortcomings indicate that the tentative adoption of the 
(leading-order) perturbative QCD Wilson coefficients for the analysis of the holographic power corrections ((3T|) is 
questionable. It will be revised on more physical grounds below. (Recall, furthermore, that the scalar background 
field of the soft wall (P5)) does not correspond to a A = 4 operator.) 

In addition to the OPE- type power corrections of Eq. (|A2[) , the holographic soft-wall correlator ([3T|) contains a two- 
dimensional power correction (times a logarithm) which cannot appear in the OPE since QCD lacks a corresponding 
(gauge-invariant and local) composite operator. However, a two-dimensional power correction of exactly this type 
was advocated some time ago and argued to improve QCD sum rule results in several hadron channels [43 |. More 
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specifically, when (possibly renormalon-related) linear contributions to the heavy-quark potential at short distances 
are approximately accounted for by an effective gluon mass A, the latter produces the correction (43 | 



ft'""'UQ^)=->^ln^(^l + 6^ + ...j (37) 

to the leading logarithm of the glueball correlator which has precisely the form of the second term in Eq. (pij) . The 
appearance of this term supports previous arguments which tentatively relate the quadratic behavior of the soft-wall 
dilaton background field ([^5]) or alternatively of A (z) [3, [3, [H, HB] to a two-dimensional power correction and 
possibly to a two-dimensional nonlocal gluon condensate . Comparison of the A correction in Eq. (f37|) with its 
counterpart in Eq. (PT|) leads to the holographic estimate 

a' ^ ^A^ (38) 

and with the approximate identification A ~ \/2Aqcd further to A^ ~ 0.15 GeV^ which is indeed of the expected 
magnitude ^4^. However, as in the case of the leading OPE power corrections the sign turns out to be opposite to 
QCD expectations, i.e. the square mass p8p is not tachyonic. 

The complete reproduction of the dependence contained in the QCD short-distance expansion (to leading 
order in as) by the soft- wall dynamics, albeit with the signs of at least the leading power corrections opposite to 
QCD expectations, suggests an interpretation which may help to disentangle the holographic predictions for Wilson 
coefficients and condensates even though they appear as products in the power corrections. Indeed, the dimensions of 
the condensates are generated by the operators of the OPE which in turn are renormalized at relatively small scales 
/i ~ 1 GeV and hence IR dominated. This makes it likely that the condensate part of the OPE and consequently 
the form of the power corrections and their scaling behavior are better reproduced by the strong-coupling dynamics 
of the soft-wall model, and that the deviations from the OPE should reside mainly in the Wilson coefficients (cf. 
Sec. in Bp . The lack of perturbative Q^-dependence due to radiative O (a^) corrections (cf. App. [X| in the soft-wall 
correlator provides additional support for this interpretation. It could be further tested by extending our comparison 
of holographic correlators with the OPE to other hadron channels. Indeed, since the condensates are universal 
(i.e. channel independent) while the Wilson coefficients are not, one would expect inconsistent soft-wall condensate 
predictions in different hadron correlator channels when relying on the questionable assumption that the soft wall 
dynamics approximates their Wilson coefficients. 

Tentatively assuming that the soft-wall dynamics approximately reproduces the values of the QCD (or Yang- 
Mills) condensates, on the other hand, one may obtain holographic estimates for the Wilson coefficients. The soft- 
wall prediction for the (leading-order) perturbative gluon condensate coefficient c|g^^^''°^ = Bq, e.g., becomes with 
(G2) ~ (20/3) A4^D m and A ~ V2Aqcd, Aqcd ^ 0.33 GeV i 

^(sw) _ 8 _ 2 (QCD,lo) ,OQ^ 

^(G2) - - :;^'-'{G^) ■ i-^-^i 

This prediction is of smaller size than the QCD value and has the opposite sign. As discussed above, it is suggestive to 
attribute at least part of these discrepancies to the strongly-coupled UV regime of the soft-wall model, although the 
estimate p9p is prone to additional error sources including the current uncertainties in the QCD value of the gluon 
condensate and its sensitivity to the presence of light quark flavors. The uncertainties in the analogous predictions 
for the Wilson coefficients of higher-dimensional operators would be further increased by the less reliably known QCD 
values of the corresponding condensates. One should note, finally, that the above approximate separation of hard and 
soft (i.e. k ^ ^) contributions to the holographic predictions would not work for the gluon mass term since both the 
mass A and its coefficient receive UV contributions. 

The soft-wall correlator in its subtracted form (P^)) fails to satisfy the low-energy theorem (jA7p : Eq. ([55)1 (if taken 
literally) implies a finite RHS while Eq. (f29|) gives H (0) = (even before discarding the contact terms), i.e. a 
vanishing LHS. Of course this comparison should be considered naive since contact terms are renormalization scheme 
dependent and devoid of intrinsic physical meaning, but other subtraction procedures, including the subtraction of 
the conformal logarithm suggested in the original LET definition [51], would lead to the same result. In fact, the 
simple soft-wall background does not correctly represent the physics of the QCD trace anomaly on which the LET 
(|A7P is based: the AdSs metric (which is dual to the energy-momentum tensor T^j, of the gauge theory (gl] on the 
flat boundary) implies {^1^}^^^^^^ = since the AdSs Weyl anomaly vanishes (G^I, and there is no scalar background 
dual to the A = 4 gluon condensate operator which appears on the RHS of the LET and in the matter anomaly 
contribution to {T^). (The soft-wall dilaton would naively correspond to a local A = 2 operator which does not exist 
in QCD but arises in (e.g. effective dual color (6^ ) theories with spontaneously broken gauge symmetry.) 
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To summarize, it is remarkable that the soft-wall background reproduces all qualitative features of the short- 
distance QCD correlator, i.e. exactly those powers and logarithms which appear in QCD, and even the hypothetical 
logarithmic corrections due to an UV gluon mass. The signs (and sizes) of both leading power corrections differ from 
those preferred in QCD, however, which we expect to be at least partly due to the failure of the strongly coupled 
UV regime to describe the perturbative QCD Wilson coefficients. Since QCD sum-rule analyses show that results for 
ground state masses and couplings (decay constants) depend sensitively on magnitude and sign of the leading power 
corrections, it is likely that the soft-wall predictions will be contaminated by this shortcoming. 

The addition of stringy corrections to the minimal bulk action ([3]) may be a promising direction for improving the 
soft-wall description in the UV. Indeed, first attempts to allow for such higher-dimensional operators in the action 
of holographic models [H, [H, HO] show that they can generate substantial contributions to the power corrections. 
Similar operators of stringy origin, including e.g. tachyon fields or a' corrections analogous to those considered in the 
vector meson sector [s^ , can therefore be expected to improve the soft-wall prediction for the short-distance correlator 
in the scalar glueball channel. 

The comparison of the above results with those for the hard-wall correlator in Sec. IIII Al shows that the whole 
nonperturbative momentum dependence of the known lOPE (up to radiative corrections) is reproduced by the holo- 
graphic hard- and soft-wall correlators in a fully complementary fashion: while the soft-wall correlator contains all 
OPE power corrections of the types induced either by gluon condensates or by an effective UV gluon mass, the non- 
pertubative physics in the hard-wall correlator is exponential and includes a term which reproduces the behavior of 
the leading instanton contributions. This complementarity of the nonperturbative physics represented by both dual 
backgrounds is likely to persist in other hadron correlators as well (at least at distances smaller than the inverse QCD 
scale) and can be exploited for diagnostic purposes, e.g. by tracing the impact of different parts of the gauge dynamics 
on hadron observables (see below). 



IV. GLUEBALL DECAY CONSTANTS 



In the following section we obtain quantitative holographic predictions for the glueball decay constants in both 
hard- wall and dilaton soft- wall backgrounds and discuss the underlying physics. 



A. Hard wall IR brane 



The values of the glueball decay constants in the hard- wall approximation may serve as a benchmark for the results 
of more elaborate holographic duals. We calculate them directly from the normalizable solutions 0, [2^ 

ip„ (z) = Nn {rrinzf J2 {rrinz) (40) 

(where n = 1,2, 3, ...) of the massless field equation ([7]) in the AdSs slice which we require to satisfy (in addition 
to the AdS/CFT boundary condition V'n {z) ^ a.t z — e ^ 0) either Dirichlet (D) or Neumann (N) boundary 
conditions on the IR brane. The normalization constants Nn are determined by the inner product of the eigenmodes, 
i.e. by requiring J^"^ dz (R/z)^ ip'^ = 1- For Dirichlet boundary conditions ipn (zm) = one then obtains the masses 
[23 , [23j and normalizations 

"4?' = = — (41) 

while the alternative Neumann boundary conditions ip'n i^m) — yield the spectrum ^24j and normalization constants 

^(N) = ihlL^ ^(N) = . (42) 

ml^^'i?3/2z„|Jo(j-i^„)| 

Here jm,n denotes the n-th zero of the m-th Bessel function ^45*1. Although the normalization constants do not affect 
the mass spectra, they provide a crucial overall scale for the decay constants. 

From the general expression for the decay constants and the hard-wall eigenmodes (HO]) one then finds 



= 1™?, 2 — 3— = (43) 
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or more specifically for the above two IR boundary conditions 

1 R3/2 1 r3/2 
f(D) = I fl— f(N) ^ I fl_ ^44) 

V2|Jl02,„)| '^^m' " V2|Jo(jl,„)| «^m' 

The expression for /l^-* can alternatively be obtained by comparing the spectral density p9p of the Neumann hard- wall 
correlator to the general spectral representation (fT2|) . This provides a useful cross-check on our calculations. 

After fixing the overall normalization factor R?'/'^ / k by comparison with the QCD gluon loop contribution according 
to Eq. (pO)) . both masses and decay constants are given (by Eqs. ((4T|) . (|42|) and (|44|) ) in terms of only one adjustable 
parameter, i.e. the IR scale ^ Aqcd of the hard-wall geometry which has to be determined from independent 
input. The resulting quantitative predictions for /„ will be discussed in Sec. IIV CI 



B. Dilaton-induced soft wall 



In the AdSs - dilaton background (I25p . the solutions of the scalar field equation ([7]) turn into Kummer's confluent 
hypergeometric functions [25| . The spectrum-generating normalizable modes then form the subset of Kummer func- 
tions whose power series expansion truncates to a finite polynomial which turns out to be of generalized Laguerre 
type [ill, i.e. 

^„ (z) = N^X^zSF, (-n, 3, z^X^) = AT^A^.^^Lf^) {X^z^) (45) 

where n = 0, 1, 2, (a)^^ = a (a -f 1) (a -f 2) ... {a + n — 1) and iFi is a confluent hypergeometric fimction [45]. The 
ensuing restriction to discrete eigenvalues = yields the glueball mass spectrum [23] 

ml = 4 (n + 2) (46) 

and relates the mass gap mo = 2\/2A to the dilaton background scale. In contrast to its hard-wall counterparts (l4l|) 
and the soft-wall spectrum grows linearly with n and thus generates a Pomeron-type trajectory |53l. |54|. 
The normalization constants iV„ are obtained from the inner product in the eigenmode space by demanding 



R 



3 



dz\-] e-^'^'^l{z) = l (47) 



which yields 



in terms of the integrals 



Jo 



(3)„ (n + l)(n + 2)- 



(Note that iV„ oc (J„) ^^"^ — > 2 '^/^n for n 3> 3, and to a rather good approximation already for n > 3.) 

From the general expression one then obtains the glueball decay constants in the soft-wall background as 

This expression shows that the fn'"'' increase by only about 40% from n = to n = oo and approach the universal 
value /m™^ = XR^/"^ / (V2k) towards higher excitation levels rather fast, in contrast to the weak but unbounded 
increase of their hard- wall counterparts (|44p . 

Eq. ([^5]) can be checked by alternatively deriving it from the spectral density (PU)) . and the factor R^/"^ /k can 
again be estimated by Eq. ((20)) which continues to hold in the soft-wall background. The dilaton scale A will be 
approximately determined in Sec. IIV CI 
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C. Quantitative analysis 



We restrict our quantitative decay-constant estimates to the glueball ground state, i.e. to y^'™-' ^ j^hw) 
y^sw) _ f^g""\ since only fs will be of phenomenological relevance in the foreseeable future and since independent 



theoretical information on it is currently available. (The extension to higher resonances by means of formulae (|44p 
and (|49p is of course immediate.) After having fixed the correlator normalization R^/k^ according to Eq. ([20)1 in 
both backgrounds, it remains to determine the IR scale z^^^ (A) of the hard- (soft-) wall gravity dual. In order to get 
an idea of how the uncertainties involved in different scale-setting approaches affect the decay constant predictions, 
we will discuss several alternative possibilities. 

A commonly adopted strategy for fixing the IR scale is to match the holographic ground state mass to lattice results. 
Uncertainties of this method include the still rather large scale-setting ambiguity of quenched lattice predictions fol] 
and the neglected light-quark effects (including quarkonium mixing and decay channels) which may substantially 
reduce the quenched scalar glueball masses [6^. Nevertheless, the quenched masses can serve as a useful benchmark 
for scale-setting purposes, in particular because it is not clear how far quark effects are accounted for in the simple 
dual dynamics which we consider here. 

We therefore base our first estimate on a typical quenched glueball mass ms ~ 1.5 GeV [H, [13, H^l , which coincides 
with the mass of the experimental glueball candidate / (1500) and fixes the IR scale of the Dirichlet (Neumann) hard 

wall at = 0.29 GeV (z^^^^ = 0.39 GeV) and that of the soft wall at A = 0.43 GeV. (Note that the values for 

Zm and A/\/2 are indeed rather close to Aqcd, as assumed in the qualitative estimates of Sec. IIIII ) When inserted 
into Eqs. (|44p and (|49| . these scales lead to the predictions 

f(D) _ 



/^^^ = 0.77 GeV, (50) 
/^N) = 0.87 GeV (51) 

in the hard-wall geometry and to the about three times smaller value 

fg'"'^ = 0.28 GeV (52) 

in the soft-wall background. Since both of the parameters which underly these results were fixed in the glueball sector 
and in the absence of quarks (recall that the estimate ((20)) is based on the free gluon loop), the above values are 
perhaps best associated with pure Yang-Mills theory. 

Alternatively, one can determine the value of the hard IR wall cutoff in the classical hadron sector, e.g. from a fit 
to TT and p meson properties as in Refs. [1, Hsj. The typical result is ~ 0.35 GeV and yields 

fP = 0.93 GeV, (53) 
/^") = 0.78 GeV. (54) 

The corresponding ground state glueball mass predictions are then mf^ = 1.80 GeV and mf^ = 1.34 GeV (where 
ms = nil). The latter is significantly smaller than most quenched lattice results but close to the /(1270) and to 
results of K-matrix analyses of scalar resonance data [gI], mixing schemes with only one 0^"*" multiplet below 1.8 GeV 
[sgj . a topological knot model [T^I and the QCD sum rule prediction ms = 1.25± 0.2 GeV One might speculate 
that fixing z^^ in the flavored meson sector takes some light-quark effects into account and hence corresponds to 
a lower, unquenched value of the scalar glueball mass (at least under Neumann IR boundary conditions). For an 
alternative estimate of the soft- wall IR mass scale A (which has not yet been determined in the meson sector), finally, 
one can use its approximate relation A ~ \/2Aqcd — 0.49 GeV (cf. e.g. Ref. [l3|) to the QCD scale Aqcd ~ 0.33 
GeV 3 (for three light quark flavors). This yields the soft- wall prediction 

4™^ = 0.31 GeV (55) 

which is similar to the first soft- wall estimate ([52)1 but corresponds to a significantly smaller glueball mass m^™-* = 1.37 
GeV. 

The above results may be summarized as follows: (i) whereas the hard-wall results for the ground state mass can 
differ by more than 30% for Dirichlet vs. Neumann IR boundary conditions, the decay constant predictions remain 
in the smaller range 

~ 0.8 - 0.9 GeV, (56) 
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and (ii) the soft-wall results for the ground state decay constant center consistently around less than half of the 
hard- wall value, 

4""^^ ~ 0.3 GeV. (57) 

The substantial difference between the hard- and soft-wall predictions can be traced to the different slope of the 
normalized dual modes at the UV brane (i.e. for z = e — > 0). (An analo gou s but less pronounced difference between 
the slopes of hard- and soft-wall modes was found in the rho meson sector [131 •) The larger slope of the hard- wall mode 
translates into a larger Bethe-Salpeter amplitude at the origin and hence into a smaller size of the scalar glueball. 

In view of the sign problem which afflicts the leading nonperturbative contributions to the soft-wall glueball corre- 
lator at distances larger than the inverse QCD scale (cf. Sec. IIIIB|) . and because of the exceptional size of the missing 
exponential contributions, one would expect the soft- wall results in the spin-0 glueball sector to be less reliable than 
their hard-wall counterparts. This expectation is corroborated by the first (quenched) lattice simulation of glueball 
decay constants [s^ which finds J^'^'^ = 0.86 ± 0.18 GeV. This lattice result is inside errors fully consistent with the 
lOPE sum-rule value /^^op^) = 1.050 ± 0.1 GeV [H, the instanton-hquid modes result f^g^^^ = 0.8 GeV ^ and 
our above holographic hard- wall result ([55]) . The soft-wall result (|57p . on the other hand, is clearly incompatible with 
the lattice prediction. 

Further insight into the holographic glueball dynamics can be gained by interpreting the above results on the 
basis of the structural complementarity between the nonperturbative physics accounted for in the soft- and hard-wall 
correlators (i.e. power vs. exponential contributions, cf. Sec. IIIip . Since the large exponential contributions to the 
hard-wall correlator can at least partially be associated with small-scale instantons and are absent in the soft-wall 
correlator, one infers that the instanton contribution can more than double the value of the decay constant. The 
mentioned lOPE sum rule analyses [sl], [s^l arrived at the same conclusion. Moreover, even the perturbative and 
hard instanton contributions alone (i.e. without the unusually small power corrections and thus analoguous to the 
hard-wall physics) were found to provide reasonable approximations to the QCD 0"^^ glueball sum rule results for the 
ground state mass and decay constant [3l[ . The neglect of the hard instanton contributions, on the other hand, leads 
to the substantially smaller prediction fg~'^^^ = 0.390 ± 0.145 GeV [7l[ which is consistent with the soft-wall result 
(|57l) but not with the lattice value. 

V. SUMMARY AND CONCLUSIONS 

We have analyzed the scalar glueball dynamics contained in two approximate holographic QCD duals, viz. the 
hard-wall IR brane geometry and the dilaton soft-wall background. Our study focuses on the 0++ glueball correlation 
function and its spectral density for which we have obtained closed analytical expressions in both gravity duals. 
A systematic comparison with the QCD physics content of the instanton-improved operator product expansion, 
a dilatational low-energy theorem and an additional, two-dimensional power correction then provides several new 
insights into the holographic representation of hadron physics as well as estimates for various bulk parameters of the 
QCD vacuum and predictions for the glueball decay constants. 

In both dual backgrounds the spectral densities are found to be non-negative, in agreement with general principles, 
and to consist of an infinite sum of zero-width glueball poles, as expected in the limit of a large number of colors. In 
their representation of specific nonperturbative glueball physics (at momenta larger than the QCD scale), however, 
both holographic duals turn out to complement each other in a mutually exclusive fashion: the soft-wall correlator 
contains all known types of QCD power corrections (to leading order in the strong coupling), generated either by 
condensates or by an effective UV gluon mass, while sizeable exponential corrections as induced by small-scale in- 
stantons are found in the hard-wall correlator. (This complementarity may in fact suggest to combine brane- and 
dilaton-induced IR physics into improved QCD duals.) 

As a consequence, the soft-wall correlator provides holographic estimates for either the three lowest-dimensional 
gluon condensates or their Wilson coefficients, as well as for the effective gluon mass (potentially associated with a 
two-dimensional nonlocal "condensate"), whereas the hard- wall correlator allows for predictions of the two leading 
moments of the instanton size distribution. All holographic estimates turn out to be of the order of magnitude 
expected from QCD, which is at least partly a consequence of the fact that the IR scale of both dual backgrounds is 
set by Aqcd- The predicted signs of the two leading dilaton-induced power corrections, however, are opposite to those 
of standard QCD estimates (and in confiict with the factorization approximation for the four-gluon condensate). We 
have argued that these shortcomings provide evidence for the short-distance physics in the OPE Wilson coefficients to 
be inadequately reproduced (beyond the leading conformal logarithm) by the strongly-coupled UV regime of bottom- 
up models. In conjunction with the absence of the sizeable exponential contributions, this casts particular doubts on 
soft-wall results for glueball observables. 
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A second main objective of our analysis was to provide first holographic estimates for the decay constants of the 
0"^^ glueball and its excitations, which contain valuable size information and are of direct importance for experimental 
glueball searches. Our analysis shows that the decay constants probe aspects of the dual dynamics to which the mass 
spectrum is less sensitive, and thus provide a new testing ground for the development of improved QCD duals. The 
hard- and soft-wall predictions for the ground-state decay constant fs differ by more than a factor of two, as do 
the corresponding QCD sum-rule results with and without hard instanton contributions. In fact, as in the sum-rule 
analyses the enhancement of fs and the consequently reduced size of the scalar glueball in the hard-wall background 
can be traced to the strong instanton-induced attraction (over relatively short distances of the order of the average 
instanton size) which the exponential contributions to the hard-wall correlator generate. It is remarkable that the 
simple hard-wall approximation can reproduce these small-instanton effects, which are known to be exceptionally 
strong in the glueball correlator. Their absence and the other shortcomings mentioned above render the soft-wall 
predictions for the glueball decay constants unreliable, while the hard-wall prediction y^'^'^^ ~ 0.8 — 0.9 GeV agrees 
inside errors with lOPE sum-rule and lattice results. 

The above arguments for the instanton-induced origin of the decay constant enhancement provide an example for 
how the complementary nonperturbative physics in the hard- and soft-wall backgrounds, which should for the most 
part generalize to other hadron channels, may be exploited to trace differences in the holographic predictions of both 
backgrounds to different origins in the soft gauge dynamics. The absence of instanton contributions to the soft-wall 
correlator provides another example: since the soft-wall dilaton background was designed to reproduce the linear 
trajectories of excited mesons, it indicates that instanton effects are not directly involved in the underlying flux-tube 
formation, in agreement with QCD expectations. 

Our results demonstrate that the comparison of holographic predictions with QCD information at the correlator 
level can provide very specific and quantitative insights into the gauge dynamics which different dual backgrounds 
encode. This holds in particular for comparisons with the QCD operator product expansion. Owing to its ability 
to factorize contributions from short- and long-distance physics to gauge theory amplitudes, the OPE allows for 
a transparent analysis and systematic improvement of several typical shortcomings of holographic models, including 
those which are rooted in their strongly coupled UV sector. These limitations notwithstanding, the amount of glueball 
dynamics which we found to be represented by even the simplest holographic duals is encouraging and indicates that 
the bottom-up approach may indeed provide a viable and systematically improvable approximation to holographic 
QCD. 

Acknowledgments 

This work was supported by the Brazilian funding agency Fundagao de Amparo a Pesquisa do Estado de Sao Paulo 
(FAPESP). 



APPENDIX A: SYNOPSIS OF QCD RESULTS 



The instanton-improved operator product expansion (lOPE) 

n(ioPE) (g2^ ^ n(oPE) (0^) 



(Al) 



of the scalar glueball correlator, which holds at spacelike momenta = —q^ ^ ^qcd; currently known up to 
operators of dimension eight, radiative corrections to the Wilson coefficients up to Oja^V, and small-size (or "direct") 
instanton contributions of O [ff') to the Wilson coefficient of the unit operator [3l|, [S^]. The standard part, with 
purely perturbative coefficients, has therefore the form (cf. [3^ . Fzij and references therein) 



n(OPE)(g2 



+ 





f Ailn 




+ A2 \n' 
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(A2) 



The full set of coefficients Ai - Di can be found in Ref. 32]. Those needed for comparison with the holographic results 
below are Aq = - {N'^ - l) / {4tt^) and (for the number of colors (hght flavors) (Nf) = 3 [U) Bq = 4 + 49q!s/ (37r), 
Co = 8 (where a small anomalous dimension correction has been neglected) and Dq = Sirag. The gluon condensates 
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are defined at the OPE scale as 



{gG') := {gfabcGl,G'>;G'^''n, 



{G') ■■ =l^{{fabcGl^Gif)-{{L,,.Gl,G';,,f) 



(A3) 
(A4) 



Contributions from instantons larger than the inverse OPE scale are accounted for in the condensates. Small-scale 
(or "direct") instantons (and anti-instantons) contribute to the Wilson coefficients, on the other hand, and affect 
dominantly the coefficient of the unit operator [3l|, [12, [Z^l ■ glueball channel, the latter is given by [H, [T^I 



i+i 



(A5) 



{K2 is a McDonald function 45|) where p and njir (p) denote the size and density of small instantons with p < in 
the vacuum. The nonpcrturbative contributions (|A5p arc known to be particularly important in the spin-0 glueball 
channels, i.e. comparable to the contributions from the perturbative coefficient and of equal or larger size than the 
power terms at > ^qcd- The expression (IA5[) can be approximated as 



(A6) 



where we specialized the instanton density to the spike distribution n^ir (p) = C*^*^ (p ~ p) which becomes exact at 
large Nc and where p and n are the average instanton size and density in the vacuum. The coupling as/ir ~ 0.2 is 
fixed at a typical instanton scale and the factor ^ ~ 0.66 excludes contributions from instantons with p > p^^ 32]. 

Further QCD information on the behavior of the glueball correlator is available in the complementary limit 0. 
Indeed, the value of the correlator at zero momentum transfer is governed by the low-energy theorem (LET) ^51i] 



n(o) 



327r 
asbo 



(G2)+0(m,) 



(A7) 



where bo — llA^c/3 — 2Nf/3, irig are the light quark masses for flavor q, and UV renormalization of both sides by a 
dispersive subtraction of high-frequency field contributions is implied [51| . The appearance of the gluon condensate 
in Eq. (|A7p reflects the fact that the LET is a consequence of the anomalous Ward Identity for the QCD dilatation 
current. Additional information on the glueball correlator has been obtained from several versions of the instanton 
liquid vacuum model (ILM) in Ref. [30], whereas direct lattice information on the (point-to-point) correlator seems 
currently not to exist. 
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